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Задача Алексідзе для відновлення потенціалу сили тяжіння 
(Представлено академіком НАН України В.І. Старостенком) 


Сформульована нова нелінійна гранична задача гравіметрії, названа задачею Алексідзе. В постановці 
А.В. Чорного вона зводиться до задачі відновлення потенціалу сили тяжіння за значеннями модуля йо- 
го градієнта на границі Ляпунова. Розглянуті деякі аналітичні властивості задачі. Її розв'язок у вигляді 
потенціалу простого шару є послідовністю розвязкі зовнішніх граничних задач Неймана для рівняння 
Лапласа за умови, що розв'язок не дуже відхиляється від заданого. Густина простого шару визначаєть- 
ся з розв'язку інтегрального рівняння Фредгольма 2-го роду для сили тяжіння. 
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Однією з непересічних проблем тлумачення даних аномалій потенціальних полів є ство- 
рення надійної методики для побудови на їх основі адекватних вимогам сьогодення і геоло- 
гічно змістовних моделей глибинної будови земної кори. Крім того, вирішення окремих при- 
кладних завдань гравіметрії 1 геодезії, пов'язаних з вивченням фігури, внутрішньої будови 
Землі, чи проявів її зовнішнього гравітаційного поля, потребує відомостей про розподіл зна- 
чень потенціалу сили тяжіння чи модуля його градієнта. Для вирішення цих завдань накопи- 
чено достатньо геофізичних даних, проте недостатньою мірою розвинуті методи видобутку з 
них всієї повноти необхідної інформації. Наявні методи трансформації потенціальних полів 
П, 2| довершують теорію потенціалу і є інформативними у вивченні внутрішньої будови пла- 
нети. Однак вони дають змогу розв'язувати з достатньою точністю ту чи іншу задачу лише за 
умови задання вхідної інформації в локальних областях певної малої міри. Спроби стикувати 
результати локальних розв'язків при відновленні інформації в глобальному масштабі зазна- 
ють невдачі |3| через відсутність точних граничних даних для розв'язання відповідних гра- 
ничних задач - Діріхле, Неймана, Стокса-Молоденського - для рівняння Лапласа. Наразі не- 
має змоги й прямо вимірювати значення гравітаційного потенціалу, однак доступні дані гра- 
вімагнітних спостережень, що являють собою значення приростів модуля градієнта потенціа- 
лу сили тяжіння (МГПСТУ). Варто скористатися цими даними при розробці схем трансформа- 
цій потенціалу в глобальній області. 

Серед трансформацій потенціальних полів чільне місце посідає задача наближеного ана- 
літичного продовження потенціалу сили тяжіння: в моделі геологічного середовища набли- 
жено задано або диференціальний оператор, або граничні умови. На практиці найчастіше для 
вирішення цього завдання використовують модель з точним диференціальним оператором 
(гармонічна апроксимація значень сили тяжіння з (х)) в задачі Молоденського |4). Відно- 


вити з гарантованою точністю поле сили тяжіння у зовнішньому просторі за його гармоніч- 
ним наближенням можна шляхом розв'язання відповідної граничної задачі в області малої 
міри, оскільки у вищезгаданих граничних задачах дані спостережень слугують лише набли- 
женими граничними умовами через негармонічність відповідних функцій. Через негармоніч- 
ність оператора трансформанти |5) (принаймні, на суходолі) різко зростає похибка визначен- 
ня потенціалу 8. (У) зі збільшенням розміру локальної області. Відтак, трансформації на ос- 
нові цих задач мають гарантовану точність лише в областях малої міри (як правило, І?х І? із 
точністю до невизначеної сталої, що залежить від геометрії цієї області). Отже, гарантована 
точність розв'язку суттєво залежить від геометрії (рельєфу і розмірів) локальної області |6), а 
критерії для поєднання локальних розв'язків в глобальних побудовах відсутні. 

Крім того, в класичному способі визначення гравіаномалій 7, (х) не враховується та обс- 


тавина, що в точках земної поверхні вектори реальної і нормальної сили тяжіння можуть ма- 


ти різні напрями через розходження поверхні земного рельєфу та референц-еліпсоїда (рис. 1), 
що особливо яскраво проявляється на стику континент-океан та в гірських районах. Просто- 
рову орієнтацію гравіметрів (нахили приладів, зумовлені різним ступенем кривизни еквіпо- 
тенціальних поверхонь, що проходять через задані пункти вимірювань) характеризує приріст 
кута а, - сов(п зи: ) з сок(п. и та) між нормалями до земної та еквіпотенціальної поверхонь 


в точці вимірювань. 





Еквіпотенціальна поверхня ху) | 





Рис. 1. Розходження векторів реальної і нормальної сили тяжіння. 
Ні цю, ані жодну іншу величину, скажімо, значення напрямних косинусів 


со5(п, х,), і-«13, що визначають напрям сили тяжіння, не вимірюють через складність 


організації таких спостережень в польових умовах. Розходження! (приріст) цих напрямів в 
реальних умовах коливається від 0" до 40" і може призвести до того, що глибинні 
неоднорідності, які зумовлюють ці розходження, не відображатимуться в класичних 
аномаліях. Необхідність вивчення нелінійної граничної задачі відновлення потенціалу сили 
тяжіння за значеннями модуля його градієнта продиктована практичною непридатністю 
класичних схем відновлення потенціалу в глобальній області. 

В зв'язку з переходом до глобальних побудов густинних моделей земної кори, які вико- 
нуються на основі даних регіональних спостережень, виникла потреба в пошуку способів 
розв'язання задачі аналітичного продовження аномалій сили тяжіння, сформульованої в |З. 
В цьому напрямі в праці |6| виокремлено два альтернативні напрямки, що базуються на ана- 
лізі характеристичних властивостей модуля градієнта потенціалу сили тяжіння. 

Перший з них потребує відшукання такого диференціального оператора, що анулює зна- 
чення модуля градієнта потенціалу поза областю розташування тяжіючих мас, та розв'язання 
для цього оператора відповідної лінійної граничної задачі. Результатом аналізу цього підходу 
є постановка і розв'язання зовнішньої задачі Діріхле для лінійного диференціального рівнян- 
ня типу Гельмгольца (яке моделює значення МГПСТ в області, не зайнятій тяжіючими маса- 
ми) з невідомим змінним коефіцієнтом, що відповідає обраному нормальному потенціалу. 
Подібному ж рівнянню задовольняють і аномалії сили тяжіння (|5|. Цей спосіб ефективний 
для продовження аномалій сили тяжіння, однак у випадку продовження значень повного гра- 
дієнта потенціалу не дає бажаних наслідків. В |3| обгрунтована і реалізована подібна схема у 
вигляді послідовності розв'язків задачі Неймана для рівняння Лапласа, які визначають збу- 
рювальний потенціал. 

Другий шлях лежить в площині постановки такої нелінійної граничної задачі для рівнян- 
ня Лапласа, в крайових умовах якої безпосередньо задіяні значення сили тяжіння. Таку зада- 
чу вперше сформулював Алексідзе в роботі (31, а в |6| її переформульовано з граничними да- 





Розбіжність напрямків векторів реальної і нормальної сили тяжіння в точках складного рельєфу коливається до кількох 
десятків секунд і може , приховати" аномалії в сотні мГал |7. 


ними для класу поверхонь Ляпунова за умови, що відновлюваний потенціал не дуже відхиля- 
ється від заданого. 

Відтак, для пошуку точніших способів відновлення наближень сили тяжіння довелось 
вийти за рамки задачі Діріхле і відшукувати уточнення коефіцієнта рівняння, первісно обчис- 
леного для нормального потенціалу, що призвело до побудови послідовних наближень по- 
тенціалу за граничними значеннями модуля його градієнта. Перехід до задачі відновлення 
потенціалу усунув необхідність обчислення наступних наближень як коефіцієнтів рівняння 
сили тяжіння, так і самих значень сили тяжіння, оскільки останні тепер можна знайти не ли- 
ше з розв'язання задачі Діріхле для рівняння сили тяжіння, а і (що простіше) з безпосеред- 
нього диференціювання відновленого потенціалу. 

У зв'язку з цим задача відновлення потенціалу за значеннями МГПСТ, яка не належить 
до класичних задач гравіметрії, набуває особливої ваги в колі обернених задач теорії потен- 
ціалу. Один з можливих способів її вирішення розроблений в праці |7|, а інший, під назвою 
»гранична задача Алексідзе для рівняння Лапласа", пропонуємо вашій увазі. Уточнімо поста- 
новку цієї задачі |3). 

За предметну модель задачі відновлення потенціалу за модулем його градієнта приймемо 
просту модель Землі як абсолютно твердого тіла, близького до тіла обертання, що рухається 
рівномірно вздовж орбіти, обертаючись навколо осі з постійною кутовою швидкістю (без 


прецесії і нутації), Нехай у -- обмежена область точок тривимірного евклідового простору, 
зайнята масами Землі, у" - необмежене доповнення до у , вільне від гравітуючих об'єктів, 
ду- границя множин у |і у , що ототожнена з фізичною поверхнею Землі. В прямокутній 
декартовій системі координат Ох, х,х, з початком у центрі Землі, осі Ох,, Ох, лежать в еква- 


торіальній площині, а вісь Ох, співпадає з віссю її обертання. Точки області у  позначмо 
грецькими літерами, б - (2, раг 2)є у ,аїї доповнення у - латинськими х - (х, зоое у 


Маси всередині Землі М (2) Сєу згустиною АМ (г) с(еаЄ генерують у навколишньому 
просторі поле сили тяжіння, що матиме потенціал 
Ох), хеу" 


0, хеу" 


ДРУ КОРА (1 
У їх Ми 2 

де / - гравітаційна стала, О(х)- 0.5а) (х, з ря - потенціал центробіжної сили, ш? - модуль 
вектора кутової швидкості Землі. Напруженість поля (значення модуля градієнта потенціалу 
за |41) дорівнює 

5, СМ(х) 0. Йо 

«-нУЖОЇ- У Яв)дхіо, о 

ув дх, дп 

де дк) ов, Хі , К «1,43 - напрямні косинуси одиничного вектора л(х) внутрішньої 


дп 


нормалі до еквіпотенціальної поверхні «ПУ (х.): (у) з Сх, яка проходить через точку х. 
Конкретизуймо постановку задачі: необхідно знайти функцію МУ), х є у", яка задо- 
вільняє всередині необмеженої замкнутої області у - у" /ду рівнянню Лапласа 


АУУ(х) «0, х є у", якщо в будь-якій точці ляпуновської границі бу області і в нескінченно 


б . 3, ММ) п. 
віддаленій точці вона задовільняє умовам: 2 й (х), хєду, У(х)-з 0 при Їх -»осо, 
ст бХ, 
де «(х) - задана неперервна функція. Називатимемо надалі цю нелінійну граничну задачу 
ззадачею Алексідзе для рівняння Лапласа". 
Проаналізуймо аналітичні особливості функції «(х) модуля градієнта потенціалу; од- 
нією з її найважливіших характеристик є наступне твердження (| 10). 
Модудь градієнта потенціалу сили тяжіння не задовольняє рівнянню Лапласа в жод- 


ній точці області у, що випливає з подання (2) і є наслідком наступної леми. 

Лема. Добуток и(х)- ц(х)- у(х) двох функцій и(х) і у(х), хе у класу с? (у) гар- 
монічний в області у тоді і лише тоді, коли кожний із співмножників гармонічний в у ,а 
їхні градієнти Уи(х) і Уи(х) ортогональні один одному в області у 

Доведення очевидно випливає з рівності 


Ух) - їі б у(х): У"и(х) З 2(Ми(х), Ум(х)) з и(х): Ух). 


ка ОХ, 
ОМ (х) 
дх, 


(о) і соб(п, рей У; Кк - 1,2,3 - ні, в чому легко переконуємося з безпосереднього обчислення. 


Кожна зі складових градієнта и(х) - є гармонічною функцією, а функції 


Якби на поверхні Землі (рівняння якої вважаємо заданим), крім значень модуля граді- 
єнта потенціалу «(х), х єду і внутрішньої нормалі т (х) до ду, вимірювали напрямок граді- 
єнта п (х) , задача відновлення потенціалу сили тяжіння звелась би до визначення потенціалу 


притяжіння У(х), х є у" як розв'язку зовнішньої задачі Неймана для рівняння Лапласа 


УРУ(х) 20, хе у", з о х єду, У(х) 
т 


-»юо, (3) 











ле НО ок о'Уе» собльті)- оо Доба), 


Однак напрямок нормалі п (х) (рис.1) невідомий через виняткову складність і вартість вимі- 
рювань. Граничні дані задачі відновлення потенціалу за значеннями модуля його градієнта 
4(х) в Е- МУ() , у рамках прийнятої моделі описує вираз |6| 


подав забав 98) 20, 














Ка к кт дх, 
2 72 
або 4.) з 4.4 (3 У сок(пь х, )-4(х) а 4) 
ке ОХ, й 
через те, що о(о) - 0. Зробімо деякі припущення. 
Хз 


Введімо нормальний потенціал, що генерується фіктивними масами, що дорівнюють 
масам в у , але розташовані в певному сенсі "нормально?" в деякій іншій області у, простої 
геометрії з границею дуг, яка не дуже відхиляється від земної поверхні бу. Подаймо потен- 
ціал притяжіння у вигляді суми У(х)- С(х)-- Т(х) нормального ((х) і збурювального Т(х) 
потенціалів, завдяки чому збурювальний потенціал описує відхилення реального розподілу 
мас в у від нормального. Нехай У(х) - одинична внутрішня нормаль до поверхні 
жи (у)е са Ух) - Уу (2) - модуль градієнта нормального потенціалу. Вважаємо за- 
даними напрямні косинуси со5(у, х, ), со5(х, , т) внутрішніх нормалей и(х), т(х) до повер- 

3 
хонь дУ, і ду, аразом з ними і со5(у, т)- У сов(у, х, сов(х, , т). 
Каї 

За таких припущень відновити потенціал притяжіння У(х), хє у" можна з граничної 
задачі (3) шляхом обчислення послідовних наближень ув 1.) Кк «0,1,2,....00: за знайденими з 
попереднього і-го кроку наближеннями сок(п,, х, ), Кк 1,43 напрямних косинусів нормалі 


п (х) визначаємо на границі бу за формулою (4) і -- 1-ше наближення сили тяжіння |7| 


ук 
влд-нк Зеньков 





,хеду, 





3 
сок(п,, т) У сов(п, ,х, )сов(х, ,т), хє ду, Ф., (х)- 4, (х)со5(п,, т) у(х)сов(и, т), хеду. 


Ка 
Після цього залишилось знайти розв'язок зовнішньої задачі Неймана для рівняння Лапласа 


ОТ, 


мі х) зФуа(х ) хеду, пох ) 


як потенціал простого шару мас неперервної густини ба (х ) х єду, розподілених з на повер- 





АТ л(х)-0, хе у", 





-», () 


хні ду 





Та 1, летів , ХЄ у" . (6) 


Невідому густину обчислюємо з нелінійного інтегрального рівняння Фредгольма дру- 
гого роду |бЇ: 


ЯМ цю Х, З 245. - 2Ф. а (2), хєду, (7) 


1 д 1 1 зі 
2ля дт, х-4 2 рез і -б 
Розв'язавши рівняння (7), наближено обчислимо з використанням (6) похідні потенціалу 
притягання 


ду") дих й от. (х) ОТ, (х) їх, У 
у (ні) Р - іч й іч ей м5 , 1,2,3. 8 
. 9 дх. дх. дх. дх. п х- ї біл 5 г 8) 





де к(х,є)- о ее 





, со5|: и, т) - Уа т, х,) 











/ / я / 

Із виразу (3) очевидно, що похідні наближень збурювального потенціалу визначаються 
у внутрішніх точках області у ,ана її границі бу значення похідних з (8) знайти неможли- 
во через наявність у підінтегральних функціях сильних неінтегровних особливостей. Для 
продовження обчислень слід знати значення похідних збурювального потенціалу саме на 
границі ду, і для їх обчислення слід передбачити спеціальну регуляризацію інтегралів (8). 
Припустимо, що таку операцію виконано 1 значення похідних (8) знайдені в точках хєду, 
тоді обчислимо наступні наближення 


3 2 
Зб а аа во о У аа) (9 ен знево 


Ка! 





Фо (х) "Чо (хсов(п, , т) - у(х)сов(м, т), хе ду 7 
після чого знову розв'язуємо граничну задачу (5)-(7), наростивши індекс і -- 1-го наближення 
(х) і наближення У ""(х)- (х)- Т, 


збурювального потенціалу Т, 2 (х) потенціалу притя- 


іч 
жіння 1 т.д. 

Строге обгрунтування заміни коректної задачі (5) розв'язком граничної задачі (6)-(7) і 
збіжності наближень УФО) до потенціалу притяжіння У(Х.), х є у" дано на прикладі близь- 


кої задачі |4). Однозначність її розв'язку доводиться теоремою єдиності розв'язання задачі 
Неймана для рівняння Лапласа за модулем його градієнта за допомогою потенціалу простого 


шару, що зводиться до доведення збіжності обчислюваних наближень УФ), К «- 0,1,2,....00 
функції У(х), хє у". Аїї, в свою чергу, легко довести, якщо виявити збіжність послідовності 
(т, (х)) наближень збурювального потенціалу: якщо Т, (х) збігається до Т(х) при К -»с0, зві- 
дси випливає не лише збіжність УЗ(х)-зу(2), х є у", але і збіжність до своїх границь будь- 


яких інших наближень, що однозначно визначаються за Т, (х). 


угіаї | р 
о мо як квадрат відношення модулів градієнтів 


збурювального та нормального потенціалу. Тоді справедлива наступна теорема |З. 





Надалі визначимо величину є (х) 


Теорема 1. Якщо величиною є" (х) можна знехтувати порівняно з «(х), то послідов- 
ність розв'язків (1, (х)) граничних задач (5) збігається до збурювального потенціалу Т(х) 
області у . 


Доведення. Досить довести збіжність послідовності 16, (х)), з якої, згідно (б), випли- 


ває збіжність (7, (х)). Оператор Дб, - 6. (х)-н Гк (х, 205, (2)45 : рівняння (7) за теоремою існу- 
бу 


вання розв'язку має обмежений обернений оператор А, А" «с«оо, звідки очевидною є 
5. )- б, (Со «2|Ф,. (х)- Ф, (2) , а з урахуванням визначення функцій Ф (х) маємо 
бі (9 - б, (9| | 20 Чіа (9| |сок(л, , т) б сов(т, 1 , т)| їі 2сЇсов(п, , т)| Чіа (2 794; (9| х 
Першим доданком у правій частині цієї нерівності можна знехтувати. Дійсно, оскільки 
7 У 9ц(х) ОТ, (9 і 2 9Т, (9 2 
-(х)- : же 2 аеМТ, ; 
9-3 Я ТО руно нут 


то, переписавши вираз для градієнта функції Т,, (х) з урахуванням очевидного співвідно- 

















оцінка 




















шення 


бо - 4. (х)сок(п, у, у) - У(х), отримуємо 4; (х)- До паї) за (дян (п, м). 


3 42 
Однак со(п, и у) щ Усов(п; | Ж )соз(х, у) 21-у7 про , тобто кут сокп, , у) за 
Ка! У 





умов теореми малий, і попередній вираз із точністю до є" (х) набуває вигляду 


дн у" оранка | (9) 





ду 
Оскільки з означення з урахуванням наближеної рівності (9) маємо 


со5(п, : т) а одна Я ти соз(х, и т) шщ со5(у, т) щ Ї -у? «Уечаі | й (0) 


Кеї К ди 








то, очевидно, що різниця со5(п, ,т)- сок(п, ,т) за порядком не перевершує числа «(х). На 
основі (9) і з урахуванням граничної умови задачі (5) матимемо 

Чіа (2) 74; (2) - со5(у, та: (х)сов(п, , т) "Чіа (х)сов(п, з , т)| , 
звідки, врахувавши умову теореми та рівність (10), одержуємо наступний ланцюжок нерівно- 


стей 
Й 














Чіа (х) 74 ОЇ - |соз(м, т) 4: (2) "Чі ОЇ з 
« |сов(и, т) а, ї (2)-- 4» (| 2 о 2 |сов(и, т) Ци, (х)-- мо (о. 
зад в зовн) 


іні , Яка Пе- 
реконує, що послідовність (6, (х)| швидко збігається в собі за умови |сок(у, т) «1; зі збіжно- 


























Остаточно матимемо нерівність 


сті 5 і (х)) в собі легко встановити збіжність до єдиної границі б (х). За умови |со5(у, т), «1 
(напрямок внутрішньої нормалі у(х) до еквіпотенціальної поверхні 0Ї/, збігається (або про- 


тилежний) з напрямком нормалі т (х) до поверхні Землі ду ) задача зводиться до зовнішньої 
задачі Неймана для рівняння Лапласа. 


Поширити формули (8), справедливі для внутрішніх точок області у , на граничні точ- 
ки хєдбу допоможе наступне твердження. 
Теорема 2. Якщо густина потенціалу простого шару (6) - неперервна на границі бу 
функція, граничні в: дае похідних 1-го порядку від потенціалу дорівнюють 


о РИ (ву боса ть). 





Доведення теореми випливає негайно із зображення потенціалу у вигляді 


года, Лани нд р в 


якщо диференціювання по х, тимчасово замінити диференціюванням по внутрішній нормалі 





т, до поверхні бу у точці х, і врахувати відому тотожність 





і, хеу 
З - - 45. з 0.5, хе у" 
Для" дт, ре - 4 Й 

0, хеу 


Індекси, якими у (3) описано послідовність граничних задач, при обчисленні похідних 
опущені, як несуттєві. Вказана в теоремі 2 формула непридатна для практичного обчислення 
похідних (через складність земного рельєфу), тому замість неї варто використовувати 
еквівалентну формулу 


ЙО ядра - Е- 48, бок т). 


Зауваження: останнє рівняння, як і попередні, що фігурують в теоремах | 12, нелінійні; 
однак, якщо зафіксувати геометрію контактної поверхні дУ(), що задана на класі Ляпунова 





С ог), то задача знаходження густини потенціалу простого шару стає лінійною і матиме 
однозначне розв'язання. 

Відтак, сформульована нова нелінійна гранична задача Алексідзе для відновлення по- 
тенціалу за значеннями модуля його градієнта і вказана послідовність її розв'язку у вигляді 
потенціалу простого шару (6), густину якого відшукують з рівняння (7). Надалі потрібно об- 
грунтувати коректність постановки цієї задачі з граничними даними на ляпуновському класі 
со) (у а . 

Автор висловлює подяку академіку В.І. Старостенку за плідні зауваження й поради, що сприя- 
ли поліпшенню рукопису статті. 
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